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Abstract. We measure the spatial distribution and hard X-ray luminosity function of cataclysmic variables (CVs) 
using the INTEGRAL aU-sky survey in the 17-60 keV energy band. The vast majority of the INTEGRAL detected 
CVs are intermediate polars with luminosities in the range 

1032-io34 erg s"^ The scale height of the Galactic 
disk population of CVs is found to be ISO^so PC. The CV luminosity function measured with INTEGRAL in hard 
X-rays is compatible with that previously determined at lower energies (3-20 keV) using a largely independent 
sample of sources detected by RXTE (located at |&| > 10° as opposed to the INTEGRAL sample, strongly 
concentrated to the Galactic plane). The cumulative 17-60 keV luminosity density of CVs per unit stellar mass 



is found to be (1.3 ± 0.3) x 10^^ erg 



Me' 



and is thus comparable to that of low-mass X-ray binaries in 



this energy band. Therefore, faint but numerous CVs are expected to provide an important contribution to the 
cumulative hard X-ray emission of galaxies. 
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1. Introduction 

The sky is much less studied in hard X-rays (>20-50 keV) 
than at lower energies because of difficulties with detect- 
ing hard X-ray photons and, more importantly, due to 
the rapidly decreasing photon flux from cosmic sources 
with increasing energy. The largest existin g catalogs com- 
prise only about ~400 hard X-ray sour ces fKriv onos et al 



most of the Galactic hard X-ray emission at energies 
below ~ 100 keV into a n umber of luminous discrete 
sources (jLebrun et al.l 120041 ). Furthemore, the remaning 



; l2007bt [bM et al.ll2007t iTueller et al.ir2007 ). As a result 
, basic properties of different populations of hard X-ray 
sources, such as number density, luminosity function, spa- 
tial distribution, etc., are poorly known. Therefore, any 
sample of hard X-ray source s with well-defined sta tisti- 
cal properties (such as that of Krivonos et al. 2007bl ) is of 
high value, since it allows one to study the parameters of 
source populations. 

The scarceness of information about the Galactic pop- 
ulations of hard X-ray sources made it difficult to inter- 
prete the results of some hard X-ray observations. For 
example, a number of early experiments detected unre- 
solved hard X-ray emission from the Galaxy and in par- 
ticular from its central region, which was suggested to 
be truly diffuse emission arising from interaction of low- 
energy cosmicj;ai[S_w^ 

Steckeilll977tlMandrou et al.|ll980l:ISacher k Schoenfeldeil 



unresolved emission - the so-called Galactic ridge X-ray 
emission (GRXE) - has been shown to closely trace the 
stellar mass distribution in the Galaxy, suggesting that it 
might be co mposed of a large num ber of fainter stellar- 
type sources ([Krivonos et al.ll2007a[ ). 



For the stellar origin of the GXRE to be correct, the 
space density of faint Galactic hard X-ray sources must be 
sufficiently high to provide the required luminosity den- 
sity. At energies above several keV, the main contribu- 
tion to the GRXE is expected from cataclysmic variables 
(GVs ), and more specifi c ally from intermediate polars 



peciii c 

qPs. ISazonov et al.lf200l iRevnivtsev et al.ll2006l ). Using 



1984 ISkibo fc Ramatv" lips'). However, the INTEGRAL 
observatory (Winkler et al. ,2003,) has recently resolved 



the lu minosity function of G Vs measured in the 3-20 keV 
band ( Sazonov et al. [looi) and assuming some fiducial 
sha pe of the spectrum of a ge neric intermed i ate po- 
lar, Revnivtsev et al. ( 2006f ) and Krivonos et al. ( 2007a ) 
demonstrated that such sources can indeed produce all of 
the unresolved Galactic hard X-ray emission. Nonetheless, 
although the spectra of intermediate po lars measured in- 
dividu all y with INTEGRAL (see e.g. 



Revnivtsev et al 



Send offprint requests to: mikej@mpa-garching.mpg.de 



l2004al lbl: iBarlow et a ll l2006l) in gener al agre e with the 
spectrum assumed by iKrivonos et al. (|2007al) . it is still 
important to directly measure the cumulative luminosity 
density of CVs at energies higher than 20 keV. 



2 



M. Revnivtsev et al.: Properties of the Galactic population of cataclysmic variables in hard X-rays 



In order to obtain direct estimates of the luminos- 
ity function and cumulative luminosity density of faint 
Galactic hard X-ray sources, we have analized results o f 



the INTEGRAL all-sky survey (jKrivonos et al.M2007bf l. 
Apart from low and high mass X-ray binaries, the 
INTEGRAL catalog contains Galactic sources of other 
types, mostly cataclysmic variables. The only non-CV 
faint Galactic sources in the INTEGRAL catalog are the 
symbiotic stars RT Cru and IGR J10109-5746/CD-57 
3057 and the peculiar stars Gamma Gas and Eta Garinae. 
Because of the very small number of such sources, we 
do not consider them here and focus on studying the 
properies of CVs detected in hard X-rays. 

2. The sample of sources 

We selected all INTEGRAL sources from the catalog of 
Krivonos etal\ (|2007bh identified as cataclysmic variables. 
This sample ideally suits our puproses because it is based 
on an all-sky survey and serendipitous with respect to GVs 
(no deep INTEGRAL exposures were devoted to observa- 
tions of CVs from the sample). We selected only those 
sources detected with more than 5a significance on the 
all-time averaged map of the sky. The resulting sample, 
consisting of 15 intermediate polars and 1 dwarf nova, is 
presented in Table [TJ 



2.1. Source distances 

Accurate distance estimates are not available for more 
than half of the sources in our sample, which presents 
a serious obstacle for constructing a GV luminosity func- 
tion. 

Usually, if the spectral type of the secondary star and 
the apparent brightness of a binary sy stem are known, 
one ca.n try to estimate i ts distance (e. g. iBarnes fc Evand 
1976t iBeuermannI 120061: lKniggel[2006h . However, in the 
case of intermediate polars (which constitute the ma- 
jority of our sample), the accuracy of this method is 
significantly affected by the fact that the optical light 
of an IP is produced not only by the secondary star 



<^nigge 200d). but also by the accretion disk 



(see e.g 
(see e.g 

face (e.g. iBonnct-Bi daud et al 



Beuermann et aD 2004) an d white dwarf sur- 
^ 200 ll ). Detailed spectro- 



scopic information allows one in some cases to determine 
the contribution of the light of the comp anion star to 
the o ptical brightness of the binary (e.g. Watson et al.l 
1995t iBonnet-Bidaud et all l200ll : iGansicke et all l2005l ). 
but such information is not available for all GVs in our 
sample. 

In order to demonstrate the importance of the contri- 
bution of the accretion disk to the NIR emission of our 
IPs, we selected those sources with available distance es- 
timates (Table [^TT|) . In all these 6 systems, the companion 
star belongs to the main sequence (we excluded here the 
IP with a giant companion GK Per and dwarf nova SS Gyg 
due to unclear correlation of its measured NIR brightness 
with the mass tranfer rate in the binary system, set by 



r d" 




1000 10" 

Luminosity (17 — 60 keV), 10^° erg/sec 



Fig. 1. Upper panel - dependence of the absolute K- 
band brightness of IPs on their orbital period. Gircles are 
measurements for sources from our sample. The dashed 
line shows the theoretical estimate of Mk based on 



Baraffe et al.l (|1998[ ). The dotted line is the same estimate 
shifted by 1.75 mag, which is equivalent to assuming that 
the secondary star contributes only 1/5 to the total near- 
infrared brightness of an IP. Middle panel ~ absolute K- 
band brightness of IPs as a function of their hard X-ray 
luminosity. An empirical fit is shown by the solid line. 
Lower panel - correlation of the hard X-ray (17-60 keV) 
luminosity of IPs with their orbital period. The dashed 
line shows the empirical dependence Lx cx P^'^. 

the orbital period) and therefore their absolute brightness 
can be related to their size. 

The orbital period (P) of such systems can be used to 
estimate the radius of the secondary star {R2) under the 
assumption that the star fills its Roche lobe: 

R2/RQ « 0.234(M2/MQ)i/3p2/3^ 

Assuming the mass-r adius relation P « 0.91 M° '^5 for the 
late-type main sequerr ce stars in CVs (e.g. IWarner 1995 : 
ISmith fc Dhillonll 19981 ). we can write ^2/^?© « O.OSlP^ i^ 
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Table 1. CVs from the INTEGRAL all-sky survey ( Krivonos et al. 2007bf ). Distance estimates enclosed in parentheses 
were estimated either assuming an absolute if-band brightness of the binary system or using the period-luminosity 
relation. 



Name 


I" 




Type 


Orb.per.(h) Fir- 


60 koV 


D,pc 




log Lhx 








V709 Gas 


120.04 


-3.45 


IP 


5.34 


3.91 


230" 


12.51 


32.54 


1.09 X 


10" 


b 


GK Per 


150.95 


-10.10 


IP 


47.9 


1.82 


470^ 


10.06 


32.80 


5.08 X 


10" 


9 


V1062 Tau 


178.08 


-10.31 


IP 


9.90 


3.67 


1100^ 


13.24 


33.87 


2.86 X 


10" 


10 


IGR J14536-5522 


319.74 


3.46 


IP 




1.20 


(347)* 
(496)* 
690^ 


12.70 


32.38 


1.74 X 


10" 


8 


IGR J15094-6649 


315.92 


-7.49 


IP 




1.05 


13.48 


32.64 


8.20 X 


10" 


9 


NY Lup 


332.44 


7.02 


IP 


9.86 


4.10 


12.52 


33.51 


7.08 X 


10" 


10 


IGR J16167-4957 


349.70 


7.33 


IP 




1.45 


(546)* 
(552)* 


13.76 


32.86 


4.27 X 


10" 


9 


IGR J16500-3307 


359.87 


8.74 


IP 




1.12 


13.71 


32.76 


5.75 X 


10" 


9 


V2400 Oph 


359.86 


8.74 


IP 


3.42 


2.64 


(331)® 
(217)* 


13.02 


32.55 


7.10 X 


10" 


9 


IGR J17195-4100 


346.98 


-2.11 


IP 




1.92 


11.69 


32.18 


3.26 X 


10" 


8 


IGR J17303-0601 


17.90 


15.01 


IP 


15.4 


3.54 


3300'^ 


12.48 


34.81 


3.47 X 


10" 


11 


4U 1849-31 


4.96 


-14.35 


IP 


3.36 


6.40 


510* 


12.64 


33.44 


8.49 X 


10" 


10 


RX J1940. 2-1025 


28.98 


-15.50 


IP 


3.36 


2.32 


230^ 


13.69 


32.31 


2.16 X 


10" 


8 


IGR J21237-H4218 


87.11 


-5.68 


IP 


7.48 


1.03 


(1854)^^ 


13.73 


33.58 


3.65 X 


10" 


10 


IGR J21335-H5105 


94.36 


-0.40 


IP 


7.19 


3.15 


(995)^ 


13.45 


33.88 


4.20 X 


10" 


10 


SS Cyg 


90.56 


-7.11 


DN 


6.60 


2.89 


166"° 


8.29 


32.12 


3.87 X 


10" 


8 


("iVlBonnet-Bidaud et al.l 


(120011), (2) 


-ISlavin. O'Brien. & Dunloi 


J (11995 


), (3)- Patterson! (| 199411. f4V 


distance 


was 


cak 



assuming K-han d absolute mag nitude Mk 
correlation, ( 7) - [Sa zonov et al. (200d), ■ 
|l995. ), (10) - iBitner. Robinson, fc B ehil 



5, (5) 



de Martino et al 



200d), based on iGansicke et al.l ^20051 ). 
fc Behil (|2007| ) 



(8) 



, (6) based on period-har d X- ray luminosity 
- iBeuermann et all (|2004l \ (;9)- l'Watson et al.l 



Table 2. Sources from the main sample with available information on distances. 



Name 


D,pc 


Orb. period, h 




log Lhx 


rriK 


Mk 








mGrab 








V709 Gas 


230" 


5.34 


3.91 


32.54 


12.51 


5.70 


V1062 Tau 


1100^ 


9.90 


3.67 


33.87 


13.24 


3.03 


NY Lup 


690^ 


9.86 


4.10 


33.51 


12.53 


3.33 


IGR J17303-0601 


3300* 


15.4 


3.54 


34.81 


12.49 


-0.11 


4U 1849-31 


510^ 


3.36 


6.40 


33.45 


12.64 


4.10 


RX J1940.2-1025 


230'' 


3.36 


2.32 


32.31 


13.70 


6.89 



( 1) -iBonnet-Bidaud et ahl tOOl\). (2) 



on IGansicke et al.M 200,5h. (5) 



PattersonI (Il994l ). (3) 
Beuermann et all (|2004l l. (6) - 



de Martino et all (|2006| ). (4) -[SazonweEaD (|2006l ). based 
IWatson et all (|l995l ) 



([Smith fc DhillonI |1998( ) , where the orbital period P is 
measured in hour s. In Fig. [1] we com pare the expected 
dependence (using iBaraffe et al.lll998l ) of the absolute K- 
band brightness of the secondary star on binary system 
period with the measured values for our CVs with known 
distances. We see that while for some sources the observed 
and predicted K-hand absolute magnitudes agree within 
one magnitude, other sources are significantly brighter in 
the near-infrared than expected. 

This deviation results from the unaccounted contribu- 
tion of emission from the accretion disk in the binary sys- 
tem. Indeed, our sources demonstrate a good correlation 
of absolute K-hand brightness with hard X-ray luminosity 
(see Fig. [1] middle panel). The linear fit to the observed 
correlation gives: 



where Lhx,30 is the source hard X-ray (17-60 keV) lumi- 
nosit}0 in units of lO'^" erg s~^ . 

This correlation is partly caused by the direct con- 
tribution of accretion disk emission to the optical 
and near-infrared light of the binary. Another con- 
tributor to this correlation is the anticipated depen- 
dence of the mass accretion rate in IPs on the or- 
bital period due to enhanced magnetic stellar wind brak- 
ing in binaries with larger optical s t ars (ISkumanich 



19721: iRappaport. V erbunt. fc Joss '19831: iPattersonI 11984 
Postnov fc Kuranov J0Q5.)- Therefore, one can signifi- 



cantly underestimate the absolute near-infrared brightness 
and distance of an IP if no allownace is made for the X-ray 
luminosity. 



Mk « 12.81 - 2.63 log Lhx, 



30) 



" Due to the relatively narrow range of spectral 
shapes of IPs in the hard X-ray en ergy band (see e.g. 
ISuleimanov. Revnivtsev. fc Ritteil l2005f ) , the hard X-ray lu- 
minosity might be considered a proxy of the source bolometric 
luminosity or mass accretion rate. 
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On the other hand, we can use the P-Lhx correlation 
calibrated on our 6 CVs with known distances to pre- 
dict hard X-ray luminosities and distances for other CVs 
from the INTEGRAL sample. The power-law slope of the 
empirical fit to the observed luminosity (mass-accretion 
rate)-orbital period correlation (Fig. [1] lower panel), 

ihx,30 ~ 9.48F3-2^ 



agrees with previous empiric al (see e.g. Patterson 



and theoretical estimates (e.g. iPostnov fc Kuranov 



19841) 



20051) . 



indicating that this dependence is physically motivated. 

We used the above empirical correlation to estimate 
the distances to those sources in our sample with mea- 
sured orbital periods. For the remaning systems we made 
a rough estimate of their distance assuming Mk = 5. 
Although this assumption is certainly inaccurate (though 
reasonable), it should not dramatically affect our analysis. 
For example, if instead of Mk = 5 we assumed Mk = 3 
or Mk = 7, the resulting value of the cumulative lumi- 
nosity density of CVs would change by less than a factor 
of 1.5. Nevertheless, one should keep in mind that since 
the distances are not known with adequate accuracy for 
all sources in the sample, all the properties of the Galactic 
CV population estimated below contain systematic uncer- 
tainties in addition to pure statistical ones. 

3. Spatial distribution 

Cataclysmic variables are not uniformly distributed 
around the Sun. Being relatively old stellar systems, they 
are concentrated towards the Galactic plane: pcv oa 
exp{—z/h), where z is the height above the Galactic plane. 
The typical scale height h of the distribution of late-type 
stars in the Galaxy is 120-150 pc. Exis ting esti mates for 



CVs fall into this range (e.g. [Patterso n 1984; lAk et al 



2008[ ) . We do not consider here the Galactocentric depen- 
dence of the CV distribution, since it only becomes im- 
portant at distances larger than 1.5-2 kpc from the Sun 
and all of our sources except for one are located closer (we 
exclude this source, IGR J17303— 0601, from the following 
determination of the scale height of the CV spatial density 
distribution) . 

Using our INTEGRAL sample of sources derived from 
an all-sky survey we can make our own estimate of the 
scale height of the Galactic disk population of CVs. To 
this end, we constructed a maximum likelihood estima- 
tor of the scale height h. The probability for a source to 
lie at distance z from the Galactic plane depends on h: 
P{z)dz = (l/h) exp{—z/h)dz. For each CV from the sam- 
ple, we calculated the maximum distance D(^l) in a given 
direction $7 to which the source is detectable by the sur- 
vey. The shape of the resulting volume V sampled by the 
survey over the whole sky was then taken into account 
in calculating the probability of source i to lie at height 
z: py {z)dz. The sum of the logarithms of such proba- 
bililities for all the sources in the sample (except for IGR 
J17303— 0601) defines our maximum likelihood estimator: 
ML = -Elni^^. 




100 200 300 400 500 

Height, pc; 

Fig. 2. Value of the maximum likelihood as a function of 
an assumed exponential scale height h of the CV spatial 
density distribution. The 68% confidence interval is indi- 
cated be the dashed lines. 

The dependence of the ML value on h is shown in 
Fig. [21 The best-fit value of the scale height is 130 pc, with 
the 68% confidence interval of 84-223 pc. This a grees with 
previous e stimates of the CV scale height (e.g. [Patterson 
1984 [Ak et al.. .2008i) . suggesting that most of our source 
distance estimates are reasonable. 

Since the h value estimated above has a relatively large 
uncertainty, we adopt h = 150 pc for estimating the CV 
luminosity function below. 

4. Luminosity function 

We employed the 1/V^max method (" Schmidt 1968[ ) to con- 
struct the hard X-ray luminosity function of CVs detected 
by the INTEGRAL all-sky survey 

To take into account the inhomogeneous distribution of 
sources around the Sun, we weight the standard (JVInax vol- 
ume found for each small solid angle 50, (at Galactic lati- 
tude h) of the survey by the space density of sources inte- 
grated over 6^1 and over distance from to dmax, the max- 
imum distance at which a given source can be detected. 
This, for an exponential distribution of sources, ppy (X 
exp(— z//t), leads to a generic voluine (iTinney et al. 1993 : 
Schwope et al.[[2002l : ISazonov et al.[l2006h 



(5K 



gon 



sm 



\b\ 



[2 - if + 2£, + 2)e 



where ^ = dmax sin Each sampled source thus con- 

tributes 1/ X] ^^cn to the estimated space density and 
l/(X]^^gcn)^ to the associated variance, where the sum 
is taken over the total solid angle of the survey. As noted 
above, we adopt the scale height h = 150 pc in our anal- 
ysis. 

The luminosity function of CVs detected by 
INTEGRAL is presented in Fig. [3] It was normal- 
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10-' 



AAul I I I I . . 

10^' 10"" 10^= 10" 10 

Luminosity (17-60 keV), erg/sec 

Fig. 3. Luminosity function of CVs detected by the 
INTEGRAL all-sky survey (solid crosses). Dashed crosses 
show the luminosity function constructed from those 
INTEGRAL CVs located at |6| > 5°, where the survey 
identification is almost complete. Arrows denote 95% up- 
per limits on the space density. 



b , I , 1 , I , , , I , d 

-20 20 

Golactlc latitude, deg 

Fig. 4. Comparison of the Galactic latitude distribution 
of CVs observed by INTEGRAL (histogram) with that ex- 
pected for an exponential distribution of CV space density, 
pcv cx exp(— z/150 pc), taking into account the sensitivity 
map of the survey. 



of this luminosity function with the all-sky one indicates 
that our determination is robust. 



ized to the local stell ar density = 0.04Mq/pc'^ 
( Jahreifi k Wielenlll997t) . We present 95% (~ 2a) upper 
limits on the space density of CVs in the lO'^^-lO'^^ and 
lO'^'^-lO'^^ erg s~^ luminosity intervals, where we have 
and 1 sources, respectively. 

The integrated space density of CVs with hard X-ray 
luminosity L17-60 keV > 10'^^ erg s~^ is pcv = (3.8 ± 
1.5) X 10"^ M-^. 

It is important to note here th at at present 



the IN TEGRAL all-sky survey catalog (jKrivonos et al 



2007bh has significant incompleteness due to the unknown 
nature of a number of sources: 25 persistent sources still 
do not have secure identifications. However, most of these 
sources reside very close to the Galactic plane and many 
of them are likely high-mass X-ray binaries, rather than 
CVs considered here. 

If there were significant incompleteness in our CV 
sample, it would manifest itself in the observed distri- 
bution of CVs over the sky. We therefore calculated the 
expected numbers of CVs in 5-deg wide Galactic lati- 
tude bins assuming the exponential density distribution 
Pcv oc exp(— z/150pc) and luminosity function N{> 



L) = L52 X 10-^(L 



17-60 kcV 



/lO'^'^erg s 



-1)-1.5 



(at 



Ln-&Q keV > 10"*^ erg s~^), and taking into account the 
sensitivity map of the survey. Comparison of the observed 
and expected distributions (Fig. |4|) suggests that our CV 
sample is not strongly incomplete. 

We can also make more quantitative assessment of pos- 
sible sample incompleteness on our results. In Fig. [3] we 
show the luminosity function constructed from the CVs 
detected away from the Galactic plane (|6| > 5°), where 
the INTEGRAL all-sky survey catalog is highly complete 
(2 unidentified sources vs. 12 CVs). The good agreement 



4.1. Comparison with the luminosity function of CVs 
detected by the RXTE all-sl<y survey 

The previous measurement of the spatial distribution and 
X-ray luminosity function of accreting white dwarfs was 
done using the RXTE slew survey at |6| > 10° in the 
energy band 3-20 keV ( Sazonov et al. 20061 ). We might 
expect that the harder energy band of the INTEGRAL 
survey should lead to preferential selection of the hardest 
sources among all sub-classes of CVs. Indeed, in our cur- 
rent INTEGRAL sample of sources the vast majority are 
intermediate pol ars, which are known to be the hardest of 
all CVs (see e.g. ISazonov et al.ll2006l and Fig. E]). 

The broadband (1 -100 keV) spectra of IPs are typi- 
cally so hard (see e.g. ISuleimanov. Revnivtsev. fc Rittei 
20051 ) that the luminosities emitted in the 3-20 and 17- 
60 keV energy bands are similar. The ratio of luminosi- 
ties in these bands does have some scatter but is typi- 
cally ii7_60 kcv/i3-20 kov ~0.6-0.9. Therefore, by divid- 
ing the luminosities of IPs measured by INTEGRAL in the 
17-60 keV by this coefficient we can compare the CV lumi- 
nosity function derived in this work with t hat determined 



with RXTE in the 3-20 keV energy band (jSazonov et al 
I2OO6I) . 

The result of this comparison is shown in Fig. [6l 
There is good agreement between the INTEGRAL and 
RXTE luminosity functions at luminosities ^3-20 keV > 
10'^^ erg s~^, where most of CVs are intermediate polars. 
At lower luminosities, other subclasses of CVs become 
numerous, but due to their relatively soft spectra such 
sources are practically absent in the INTEGRAL (hard 
X-ray selected) sample. 
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Energy, keV 

Fig. 5. Typical broadband spectra of different classes of 
CVs, indicating that INTEGRAL observations in the hard 
X-ray energy band (17-60 keV) are biased toward detect- 
ing the hardest CVs - intermediate polars. 

5. Hard X-ray luminosity density of CVs 

Using the same l/Vmax method we can also estimate the 
cumulative hard X-ray emissivity of CVs: 



EM 



5EM 



VE'5Knax(i^) 



We find that the cumulative 17-60 keV emissivity of CVs 
with Ln-Qo kov > 10^^ erg s~^ per unit stellar mass 
is L17-6O kov/M* = (1.3 ± 0.3) X 10^7 erg s'^ Mq\ 
Fig. [7] shows the cumulative emissivity as a function of 
the threshold luminosity. 

The obtained value of the CV cumulative emissivity is 
in remarkably good agreement with the emissivity of the 
unresolved hard X-ray emission observed by INTEGRAL 
from the Galactic plane reg ion: Ln-apkeV /J}-^* = (0.9 



1.2) X 10^"^ erg s'^ Mg^ (jKrivonos et all l2007al) . sug- 



gesting that this emission is mostly the superposition of 
numerous cataclysmic variables. 

6. Broadband hard X-ray luminosity function of 
the old stellar population in the Milky Way 

It has been known since early observations with focusing 
X-ray telescopes that the cumulative emission of bright 
(Lx > 10'^^ — 10'^^ erg s~^) discrete sources constitutes 
a si gnificant fraction of the total X-ray emission of galax- 
ies ()Trinchieri fc Fabbianol [l985 ) . This is especially true 
for energies >2-3 keV, at which emission from hot inter- 
stellar gas cannot contribite significantly even in gas rich 
galaxies because typical gas temperatures are kT ^ 1 keV. 



10-! 



io-» 



10-' 



J I I 

1032 1033 1034 

Luminosity (3-20 keV), erg/sec 



Fig. 6. Luminosity function of CVs (mostly intermediate 
polars) detected by the INTEGRAL all-sky survey (solid 
crosses), converted to the 3-20 keV energy band assum- 
ing the hardness ratio Li7_60kcv/-^3-20keV = 0.7 for IPs. 
This luminosity function is compared with the luminosity 
function of faint Galactic sources detected by the RXTE 
slew survey at |6| > 10° (dashed crosses), with dotted 
crosses sh owing the contributi on of IPs to this luminosity 



function ( Sazonov et al. I l2006h . Note that, in contrast to 



Fig. [21 we show here a rough estimate of the space density 
in the Lir-eo keV = 10^*-10^^ erg interval based on 
the single source detected by INTEGRAL in this range, 
rather than an upper limit. 



Therefore, it is often assumed that most of the hard X- 
ray emission of galaxies is provided by low- and high-mass 
X-ray binaries, while the contribution of weaker sources is 
usually neglected. 

However, as we have shown in the previous section, 
CVs are expected to produce a hard X-ray emissivity com- 
parable to that generated by brighter sources - low-mass 
X-ray binaries . The s ame was previously suggested by 



Krivonos et al.l (j2007ai ) based on a study of the Galactic 



bulge with INTEGRAL in which the luminosity of the un- 
resolved hard X-ray emission was found to be comparable 
to the cumulative luminosity of bright resolved sources. 
We therefore construct below the combined hard X-ray 
luminosity function of CVs and LMXBs. 

In Fig. [8] we present the luminosity function of CVs 
and LMXBs in the 17-60 keV energy band based on 
INTEGRAL data analysed in the present paper and in 
Revnivtsev et al. 112008). This plot demonstrates that the 



cumulative hard X-ray luminosity density of LMXBs is 
only marginally larger than that of CVs. 

This ratio of the cumulative emissivities of LMXBs and 
CVs in the hard X-ray band is significantly smaller than in 
the standard X-ray band (2-10 keV) see e.g. lSazonov et aL 
2006). The major reason for this is that the brightest 
LMXBs (ix ~ 10^^ erg s"! in the 2-10 keV band) which 
provide the majority of cumulative emissivity of LMXBs 
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Fig. 7. Cumulative hard X-ray luminosity density of CVs 
determined from the INTEGRAL all-sky survey. 

have very soft spectra and thus emit only a small frac- 
tion of their total emission in hard X-rays. At the same 
time, luminous CVs emit almost the same fraction of their 
bolometric luminosity in standard and hard X-rays. 

7. Conclusions 

We have studied a sample CVs detected in the hard X-ray 
band (17-60 keV) by the INTEGRAL all-sky survey. This 
sample is one of the largest unbiased samples of luminous 
intermediate polars existing today. 
Our findings are: 

1. The scale height of the Galactic disk distribution of 
CVs ish= 130t% pc. 

2. We have measured the luminosity function of CVs in 
the 17-60 keV energy band. The space density of CVs 
with luminosities ii7-60 keV > 10^^ erg s~^ is (3.8 ± 
1.5) X IO-^Mq^ 

3. The luminosity density of CVs with luminosities 
Li7_6o keV > 10^^ erg s^^ is iiy-eo kcv/M^ = (1.3 ± 
0.3) X lO^"^ erg s"! Mq^ This value is in remark- 
ably good agreement with the emissivity of the unre- 
solved (with INTEGRAL) hard X-ray emission in the 
Galactic plane region. 

We emphasize that all our estimates depend on the 
accuracy of the distance determination for our sources. 
Therefore, the current lack of accurate distances for a 
number of sources from the sample introduces additional 
(not included in the errors quoted above) uncertainties 



Fig. 8. Combined luminosity function of CVs (Lhx < 
1Q34.5 gj.g g-i^ £iigd circles) and LMXBs (Lhx > 10^^'^ 
erg s~^, open circles). 

in the measured parameters of the Galactic population of 
cataclysmic variables. 
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